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Abstract — This work proposes the use of automatic co-registered satellite images to obtain 
large, high frequency and highly accurate shoreline time series. Very high-resolution images 
are used to co-register Landsat and Sentinel-2 images to reach a vertical and horizontal spatial 
accuracy within 3 m. Satellite derived shorelines presented positioning errors lower than 
mission pixel’s resolution. A discussion is presented on the applicability of those shorelines 
through an application to Tordera Delta (Spain). 


Introduction 


Climate Change and human activity are likely to have an increasingly dramatic 
effect on coastal zones. The analysis of the effects on the coast requires better knowledge and 
more regular and efficient monitoring in order to manage and mitigate impacts of the 
changing coastal environment. Although in-situ measurements are highly efficient on 
capturing coastal parameters and features at a given time, the cost of continuous acquisition 
campaigns is dissuasive. Earth Observations provide wide spatial coverage over a large 
temporal scale. However, the lack of geolocation accuracy (i.e. the position of details and 
features usually visible on HR satellite imagery) challenges the conclusions of most academic 
studies. Spatial accuracy is defined for this study as the conformity in location between the 
different earth observations and a reference, in a way that an object in a set of images is 
located at the same location, while the precision is the spatial shift of the objects from what 
is actually observed in situ. The specific dynamics linked with the coastal environment imply 
additional constraints due to the different scales of the studied phenomena. Coastal processes, 
which impact the coastal morphology, are characterized by high temporal variability, ranging 
from a few hours or days (e.g. storms or artificial beach nourishments) to years or even 
decades (e.g. Metonic cycles and sea level rise). 

The reduced spatial and time scale of coastal changes makes high spatial accuracy 
necessary and this is under constant improvement (e.g. [7]). The application of the geo-location 
method, for example, greatly improves the spatial accuracy and allows features inter-location 
around a few meters accuracy, although such procedure can be time consuming. While the 
use of automatic methods to obtain geo-located images can be very effective, so far it has not 
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been tested with the focus on detecting coastal changes (e.g. [11]). On the other hand, it is 
the temporal variability that challenges the most coastal analysis. While there is an evident 
interest on obtaining long, high frequency and high accurate time series from EO products, 
building a stable statistical data set from them is a challenge. [14], for example, demonstrated 
that current available satellite data from optical sensors are suitable to represent shoreline 
changes in inter-annual to decadal time scales, but could not explain sub-annual variations. 

This study proposes the use of automatic co-registered satellite images to obtain 
long, high frequency and high accurate shoreline time-series. These Satellite Derived 
Shorelines (SDSL) are used to compose long and high frequency time series. The 
applicability of those products is demonstrated with a study case at Malgrat beach (Tordera 
Delta - Spain) in which shoreline variation is assessed in different time-scales. 


Materials and Methods 


Malgrat beach is located southward of Tordera river mouth (forming the Tordera 
Delta) in the northern coast of Cataluña (Spain — Figure la). The longshore sediment 
transport occurs predominantly from NE to SW in Cantalufia coast. Over years, the water 
capture along the river basin and the construction of port structures on the coast northwards 
of Tordera River (e.g. Port of Blanes), have led to a considerable reduction of the sediment 
input and have resulted in intense erosion of Malgrat beach ([9], [2]). In the 90’s, a 
considerable amount of sand was recovered from the shallows in front of the mouth of the 
Tordera river in order to re-nourish beaches in Cataluña. The resultant trough in the sea floor 
is now a sink to the sediment and makes redistribution of the sand along the coast difficult. 

In the past decades, this coast suffered a maximum shoreline retreat of about 120 m 
(Figure 1b and Figure 1c), with erosion rates reaching -3.8 m/year ([5]). Rapid retreat of 
about 25 m can also result from single storm events ([6]), a value that can be enhanced due 
to changes in wave conditions derived from Climate Change ([10]). These issues led to the 
need for numerous interventions, including soft stabilization and nourishment projects 
(Figure 2). Nonetheless, these initiatives did not prevent the erosion occurring, and a large 
urban area has been lost due to shoreline retreat. 

The historical shoreline variation at Malgrat beach is characterized by a sequence 
of retreats due to storm induced erosion events alternated with punctual human-induced 
progradation (Figure 2). In 2008, for example, the San Esteve Storm, desolated Cataluña 
coast and resulted in significant shoreline retreat at Tordera Delta [4]. In May 2015, the 
continuous erosion required the authorities to carry out a project to restore and protect the beach, 
with the construction of geotextile groins and the dumping of 112 000 m3 of sand (Figure 2b). 
In the following November, a storm event hit the Catalufia coast and a significant amount of 
the sand dumped on the beach in the previous months was lost. In January 2017, a sequence 
of storm events again made the shoreline retreat (Figure 2c), reaching a maximum on 19""-23"4 
of that month. The most recent event with significant magnitude occurred on January 2020, 
a weather event named the Gloria Storm, which damaged urban areas near the coast due to 
severe wave and storm surge conditions. As part of the plan to monitor the area, in situ 
measurements were taken, for example, on November 11" 2015, and May 2017. All the 
occurrences outlined above make this area an interesting site for EO application in different 
time scales, since change due to long-, mid- and short-term processes take place in it. 
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Maigrat beach (2004) Malgrat beach (2018) 


Figure 1 - Location of the study case site (a), Malgrat beach (Spain) and 
situation of the coast in 2004 (b) and 2018 (c). The yellow polygon 
indicates the coastal area lost due to erosion problems along the 14 years. 
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Figure 2 - Timeline of punctual actions/events that took to espressive shoreline 
variability of Malgrat beach. a) 2008: results from LiDAR analysis from [8] where red 
colors indicate sediment lost areas due to Sant Esteve Storm. b) 2015: nourishment 
project carried out to restore the beach. c) 2017: effects of the storm registered by 
local news. d) 2020: aerial photos taken of previous and posterior to the Gloria Storm. 
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For this study, we use a combination of several optical sensors to increase the 
number of observations over our area of interest and to obtain longer time series. Landsat 5 
& 8 and Sentinel-2 data were combined to cover a wider time range (1994-2018) and to 
increase the number of observations over a same site. Very High Resolution (VHR) imagery 
is used to co-register the Sentinel-2 and Landsat data. Those VHR are obtained from 
commercial missions. Maxar’s WorldView 2 data were used for the co-registration of HR 
images over Tordera area. Table 1 details the specification of the different missions used. 


Table 1 - EO missions and their specifications with optical sensors. 


Satellite sensor & bands pixel resolution of Revisit Years 
constellation (optical) L1 products time active 
Landsat 5 0.45+0.52 um, 0.52 + 0.60 um, 30m 
0.63 + 0.69 um, 0.77 = 0.90 um 
Landsat 8 (OLT) |0.503 + 0.676 um 15m 16 days 2013 — 
0.435 + 0.451 um, 0.452 + 0.512 um, 30m 


0.533 + 0.590 um, 0.636 + 0.673 um, 
0.851 + 0.879 um 


Sentinel-2 A & /448+546 nm, 537+583 nm, 10m 5 days 2015 — 
B /MSI 645+683 nm, 762+908 nm 
604+723 nm, 731+749 nm, 20m 
768+796 nm 
430+467 nm, 932+958 nm 60 m 
WV-2 450+800 nm 0.46 m 1 day 2009 — 
400+450nm, 450+510 nm, 2m 


510+580 nm, 585+625 nm, 
630+690nm, 705+745 nm, 
770+895 nm, 860+1040 nm 


Optical observations were used to compose large and high frequency time series of 
shorelines. For that purpose, the waterlines (the observed instance in time of the boundary 
between land and sea) need to be extracted with a very high location accuracy. However, 
from one sensor to another and from one acquisition to another, we observe a spatial shift 
(i.e. EO observations have errors and uncertainties that lead to different positioning). For 
Landsat-8 the geo-location accuracy is around 30 m, while for Sentinel-2 L1C (without 
ground control points) the accuracy is around 10 m at 20 confidence level [3]. To correct the 
displacement between images, VHR imagery with high location accuracy, are used as a 
master to co-registered observations of the different sensors. VHR were pre-processed to 
correct and remove the effect of haze vapor, then, the ASORICS (Automated and Robust 
Open-Source Image Co-Registration Software) [11] was used to apply the Fourier shift 
theorem in order to determine precise X/Y offset at a given geographical position through 
phase correlation. This process can only be applied for single band images with same spatial 
resolution and, ideally, similar pixel intensity. This criterion is not fulfilled when dealing 
with multi-temporal and multi-spectral data from different sensors, therefore, we combined 
a pure phase correlation approach with additional processing and evaluation modules. 

Once the satellite images were geo-located, waterlines were obtained and 
transformed to datum (Mean Seal Level REDMAR — Red de Mareografos de Puertos del 
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Estado) based shorelines using the beach slope and sea water level (observed sea level at the 
coast including tide, surge and wave data). SDSLs were validated for the study site using data 
from in-situ measurements. SDSLs closest to measurement days (maximum +7 days) were 
compared with in-situ shorelines at transects displaced every 15 m along the coast. The cross- 
shore distance between in-situ measured data and the SDSL was used as an error measure to 
calculate the Mean Absolute Error (MAE) and the Brier Skill Score (BSS, [12]). Perfect 
agreement gives a BSS of 1, whereas observing the baseline condition (first SDSL from 
1994) gives a score of 0. Additionally, shoreline variations were compared to results 
presented in previous studies. 

A detailed analysis of shoreline variation on Malgrat beach was carried out to show 
the potential of co-registered SDSL time series to assess coastal changes in different 
timescales. A total of 184 shorelines from 1994 to 2018 were used to assess short-, medium- 
and long-term shoreline variation. Long-term trends in shoreline changes were estimated 
using the Digital Shoreline Analysis System, DSAS [13]. Shoreline evolution was assessed 
to verify the changes in trends of Tordera Delta coast. Short term application was assessed 
through the analysis of the shoreline movement in two events: i) a nourishment event carried 
out in May 2015 and the subsequent storm erosion in November 2015, and ii) a sequence of 
storms on January 2017 (see list of events in Figure 2). 


Results 


The co-registration process, performed on all public images, allows the relocation 
and the perfect overlap from one EO image within the same tile to another. In 90 % of cases, 
the co-registered images showed horizontal and vertical shifts below 3 m in comparison to 
VHR. In other words, all images are repeated within 3m of their true location. 

In-situ measurements and SDSL from November 2015 are used for validation. The 
shorelines are presented in Figure 3, as well as the baseline used to calculate the BSS. Here, 
the baseline was taken as the earliest SDSL available (from 1994). The MAE between both 
shorelines was equal to 21 m, falling within the range of the horizontal accuracy of Landsat 
8, which is typically 30 m. BSS was equal to 0.96 which indicates good agreement between 
measured and satellite data and the ability to represent long-term variations. Shoreline 
variation rate, calculated using SDSLs, also agrees with values presented before. [1], used 
high-resolution aerial photographs and obtained an average erosion tendency of -4.68 m/year 
between 1995 and 2009 in the coast southward of Tordera river moth, a value similar to the 
one obtained here using the SDSL (-4.79 m/year) for the same period. 


Discussion 


The co-registration process results in images with high spatial accuracy, however 
the global quality depends on the number of ground points that the processor is able to use to 
adjust the images. Points selection, which is linked with weather conditions such as cloud 
cover and type of land use (points extraction is more reliable over man-made structures, over 
build-up areas than over natural environment such as forest or over large crops field) is a vital 
step in establishing the process. 
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Figure 3 - a) Position of the cross-shore transects, shorelines (measured and satellite) and 
reference line; b) cross-shore distances between shorelines and between shorelines and 
the reference. 


The increment in accuracy obtained through co-registration allows more reliable 
analysis of coastal changes in different time-scales. Here we show the applicability of this 
kind of product to both a long-term, (annual) and event scale. The net shoreline movement 
and the variation rate (calculated with a linear regression method) obtained from the long- 
term analysis are presented in Figure 4. Negative values indicate erosive processes. A retreat 
tendency is observed along the whole beach, with maximum of -6.19 m/year in the central 
area and minimum of -2.28 m/year in the northern zone. The most critical area shows a retreat 
of almost 150 m during the period of 25 years analysed (e.g. transect 33). The evolution of 
shoreline displacement and the yearly average of the variation rate are also represented in the 
Figure 4c, in terms of distance from the reference line, which is located inland (see reference 
line location in Figure 4a). Although some increment can be observed along the time series 
(e.g. years 1997, 2005 and 2015 in Figure 4c), a general retreat is clear over the 25 years, and 
is especially critical on transect 33. This critical zone corresponds to the area typically used 
by camping services (example of urban services occupying the dry beach) and which was 
completely removed due to long-term processes and due to extreme storm events. 
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Figure 4 - Long-term analysis (SDSL 1994-2018). Spatial assessment: a) net shoreline 
movement along the beach, b) variation rate along the beach (negative values indicate erosive 
process). Temporal assessment: c) average distance from the reference line located inland, 
d) average shoreline variation (negative values indicate erosive process). 


The most recent reported extreme events that affected Malgrat occurred over three 
distinct episodes, on November 2015, January 2017 and January 2020, and they contributed 
to the significant reduction of the beach width. The analysis of short-term shoreline variation 
(e.g. storms and nourishment) using satellite images is quite challenging due to the usual lack 
of short-term registers. In addition, clouds are frequently an issue obscuring optical detection 
during storms. While the problem related to clouds cannot be solved for shorelines derived 
from optical sensors, the increment on the frequency of acquisition and accuracy in recent 
satellite missions is a step forward and allows for the punctual analysis of events capable of 
moving considerable amounts of sand from the coast. Here, we analyze the shoreline 
evolution during 2015 and 2017 events individually by assessing the shorelines registered 
immediately before and immediately after those events. 

Shoreline advance due to the nourishment carried on May 2015 is represented in 
Figure 5a (SL movement 1: movement from March 2015 to June 2015). A maximum advance 
of about 60 m was observed and agrees with the nourishment project presented by that date 
(see Figure 4b). The following storm event (November 2015) that hit the shoreline causing 
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retreat was also detected (SL movement 2: movement from June 2015 to November 2015) 
leading to the displacement of the coast up to 40 m landward. In the following months, the 
retreat continued along almost the entire shoreline, moving up to 30 m landward in some 
areas (e.g. transects 2 to 17 and transect 57). 
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Figure 5 - Shorelines pre and post events from 2015 (nourishment and storm erosion - a) and 
2017 (sequence of storm erosions - b). Left panels show the net shoreline movement between 
the shorelines register. 


The sequence of storm events that occurred on January 2017, is represented in 
Figure 5b. The initial shoreline retreat observed from November 2016 to January 2017 (SL 
moment 1) reaches values of about 30 m in the critical area around transect 25. The erosion 
problem is enhanced in the following month (SL movement 2: movement from Jan/2017 to 
February 2017), in which the retreat occurs over the previously retreated shoreline. The 
shoreline advances over the camping area and reaches the sea front as observed in-situ (see 
Figure 2c). 
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Conclusions 


The analysis presented here showed the potential use of large and dense SDSL time- 
series, populated with products from multiple satellite missions to enable the analysis of long- 
term variations. The importance of applying the co-registration method, which improves the 
spatial accuracy, on the assessment of coastal variability was also discussed. This well-known 
methodology has been extensively explored in the past for many purposes, but the automatic 
application and its use to assess coastal changes is an important innovation. 

The increased frequency of acquisition in recent satellite missions allowed the 
analysis of short-term events to demonstrate significant shoreline variations (variations 
higher than the shoreline accuracy). Such advances are a step forward and improve the 
analysis of the effect of storm events and nourishment on shoreline variations. The possibility 
of using products from radar sensors (not explored here), which allow shoreline acquisition 
under cloudy conditions, may also grant important additional analysis of this kind of events, 
and will be further explored in a future work. 
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